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ABSTRACT 

Context. The emission line of [C n] at 158 fim is one of the strongest cooling lines of the interstellar medium (ISM) in galaxies. 

Aims. Disentangling the relative contributions of the different ISM phases to [C n] emission, is a major topic of the HerM33es program, a Herschel 

key project to study the ISM in the nearby spiral galaxy M33. 

Methods. Using PACS, we have mapped the emission of [Cn] 158 yum, [Oi] 63yum, and other FIR lines in a 2'x2' region of the northern spiral 
arm of M33, centered on the Hn region BCLMP 302. At the peak of Hey emission, we have observed in addition a velocity resolved [C n] spectrum 
using HIFI. We use scatterplots to compare these data with PACS 160/jm continuum maps, and with maps of CO and Hi data, at a common 
resolution of 12" or 50 pc. Maps of Ho- and 24 jim emission observed with Spitzer are used to estimate the SFR. We have created maps of the [C n] 
and [O i] 63 jim emission and detected [N n] 122 yum and [N in] 57 jim at individual positions. 

Results. The [C n] line observed with HIFI is significantly broader than that of CO, and slightly blue-shifted. In addition, there is little spatial 
correlation between [C n] observed with PACS and CO over the mapped region. There is even less spatial correlation between [C n] and the atomic 
gas traced by H i. Detailed comparison of the observed intensities towards the H n region with models of photo ionization and photon dominated 
regions, confirms that a significant fraction, 20-30%, of the observed [C n] emission stems from the ionized gas and not from the molecular cloud. 
The gas heating efficiency, using the ratio between [C n] and the TIR as a proxy, varies between 0.07 and 1.5%, with the largest variations found 
outside the H n region. 

Key words. ISM: clouds - ISM: HII regions - ISM: photon-dominated regions (PDR) - Galaxies: individual: M33 - Galaxies: ISM - Galaxies: star 
formation 
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1. Introduction 

The thermal balance and dynamics of the interstellar medium in 
galaxies, is best studied through spectroscopic observations of 
its major cooling lines: [Ci], [Cn], and [Oi] trace the transition 
regions between the atomic and molecular gas, while CO traces 
the dense molecular gas that provides the reservoir for stars to 
form. Most of the important cooling lines lie in the far-infrared 
and submillimeter regime. Therefore, it is difficult or impossi- 
ble to study them with ground-based telescopes, while previous 
space-based telescopes provided low sensitivity and coarse an- 
gular resolution. The infrared line emission is mostly optically 
thin and can be traced throughout the densest regions in galax- 
ies, allowing an unhindered view of the ISM. Herschel provides, 
for the first time, an opportunity to image the major tracers of the 
ISM at a sensitivity, spectral and spatial resolution that allows to 
study the interplay between star formation and the active ISM 
throughout our Milky Way and in nearby galaxies. 

The two fine structure lines of [Cn] at 158 /urn and [Oi] 
at 63 /im, are the strongest cooling lines of the ISM, carrying 
up to a few percent of the total energy emitted from galaxies 
in the far-infrared wavelengths. The [Cn] line lies 92 K above 
the ground s tate with a critical de nsity for collisions with H of 
3 x 10 3 cnT 3 dKaufman et al.ll999l) . While both lines are thought 



* Herschel is an ES A space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 



to trace photon-dominated regions (PDRs) at the FUV irradiated 
surfaces of molecular clouds, it has been realized early-on that a 
non-negligible fra ction of the [ C n] emission may stem from the 
ionized medium (lHeileslll994l) . Owing to its higher upper en- 
ergy level (228 K) and higher critical density of 5 x 10 5 cm -3 , 
the [O i] 63 /im line is a m ore dominant cool ant in warmer and 
denser neutral regions (e.g. lRollig et alj|2006l) . 

The [Cn] line together with the [Oi] line are diag- 
nostics to infer the physical conditions in the gas, its 
temperatures, densities, and radiation fields, by compar- 
ing the intensities and their ratios with predictions of 
PDR models (e.g.. iTielens & Hollenbachl 11985b IWolfire et al 
I l990|: iKaufman et al 1 119991: iRollig et al. 1 1200 
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1998). Previous observational studies of [Cn] emission 
from external galaxie s include the statistic al studies by 
ICrawford et aTJ d!985l) : IStacev et al.l dl991l) : iMalhotra et all 
(1200 ll) . More recently t he [Cn] emission from a few individ- 
ual g a laxies e.g., LMC dlsrael et al. Ill996 | ). M51 dNikola et al.l 
1 200 n iKramer et al.l 120051). NGC 69 46 (Madden et al.l 119931: 
bontursi et al.l l2002h . M83 dKramer et al.1 120051). NGC 1313 
dContursi et al.ll2002h. M 31 dRodriguez-Fernandez et al.ll2006l) . 



NGC 1097 dContursi et al.l 120021: iBeiraoet al.l 1201 Ol) have been 
observed. These papers explore the orig in of the fC i i ] emis sion. 
Of these studies, the study of LMC by llsrael et al. I d 19961) was 
at a resolution of 16 pc and iRodriguez-Fernandez et al. H2006b 
resolved the spiral arms of M 3 1 at spatial scales of 300 pc. 
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Fig. 1. 250/mi Herschel-SPIRE dust continuum image of 
the northern spiral arm of M33 along with the locations 
of the prominent Hn regions NGC604, NGC595, IC 142, 
BCLMP 302, and BCLMP 691. The rectangle delineates the 
area observed with PACS and HIFI, which is centered on 
BCLMP 302. The white circles indicate the positions and 70" 
beam of ISO/LWS observations done within the BCLMP 302 re- 
gion. The white ellipse delineates 2 kpc galacto centric distance. 
The two parallel lines running along the major axis of M33, mark 
the strip for which we plan to observe the [C n] and other FIR 
lines using PACS and HIFI. 



In addition to [C n] and [O i] 63/mi, there are several addi- 
tional FIR lines, the mid-7 CO transitions, and lines of single 
and double ionized N and O which provide information about 
the gas density (multiple transitions of the same ions), hardness 
of the stellar radiation field (ratio of intensities of two ionization 
states of the same species), and ionizing flux. In particular the 
[Nil] lines at 122 and 205 /mi allow us to estimate the density 
of the ionized gas, which is a key parameter to model the [C n] 
emission stemming from the ionized gas. 

Although much detailed information can be obtained by 
studying nearby (Milky Way) sites of star formation, a more 
comprehensive view is possible with a nearby moderately in- 
clined galaxy such as M33. M 33 is a nearby, moderately metal 
poor late-type spiral galaxy with no bulge or ring, classified as 
SA(s)cd. It is the 3rd largest member of the Local Group. Its 
mass, size, and average metallicity are similar to those of the 
Large Magellanic Cloud (LMC). M33 hosts some of the bright- 
est Hn complexes in the Local Group. NGC604 is the second 
brightest Hn region after 30Dorad us in the LMC. Its inclina- 
tion (2 = 56°) dRegan& Vogelll 19941) yields a small line-of-sight 
depth which allows to study individual cloud complexes not suf- 
fering from distance ambiguities and confusi on like Galactic 
obser vations do. Its close distance of 840 kpc dFreedman et al] 
1199 lb provides a spatial resolution of 50pc at 12", allowing 
us to resolve giant molecular associations in its disk with cur- 
rent single dish millimeter and far-infrared telescopes like the 
IRAM 30m telescope and Herschel. Its recent star formation 
activity, together with the absence of signs of recent merg- 
ers, makes M 33 an ideal source to study the interplay of gas, 
dust, and star formation in its disk. This is the aim of the open 
time key program "H erschel M33 extended survey" HerM33es 
dKramer et alJ 120101) . To this end we are surveying the major 
cooling lines, notably [C n], [O i], [N n], as well as the dust spec- 
tral energy distribution (SED) using all three instruments on 




RA (J2000) 

Fig. 2. Maps of 158 /mi [Cn] (in color and black contours) 
and 63 /mi [O i] (white contours) emission observed with PACS 
toward BCLMP 302. The [Cn] intensities shown in the color 
wedge are in units of erg s _I cm 2 sr _1 . The [Oi] 63 /mi map has 
been smoothed to 12" for easy comparison with [Cn]. The Ha 
peak observed with HIFI is marked with the asterisk. The white 
dots show the footprint of the PACS observations. The contour 
levels are at 10 to 100% (in steps of 10%) of peak [Cn] inten- 
sity of 1.18xl0 _4 ergcm _2 s _1 sr~'. The contour levels for [Oi] 
63 /mi emission are 30 to 100% (in steps of 10%) of the peak of 
3.0xl0~ 5 ergcirT 2 s sr~'. Both images are at a common reso- 
lution of 12". 



board of Herschel, HIFI, PACS, and SPIRE. We also use ancil- 
lary observations of Ha, H i, CO and dust continuum at 24 /mi. 
The HerM33es PACS and HIFI spectral line observations will 
focus on a number of individual regions along the major axis of 
M 33 which will initially be presented individually. 

Here, we present first spectroscopic results obtained for 
BCLMP302, one of the brightest H n regions of M33. We discuss 
a 2' x 2' (~0.5 kpc x 0.5 kpc) field in the northern spiral arm at a 
galactocentric distance of 2 kpc (Figure[T). This region lies along 
the major axis of M33, and hous e s one o f the bright H n regions, 
BCLM P 302 dBoulesteix et all 1 1974b llsrael & van der Krurd 
119741) . Using ISO/SWS. IWillner & Nelson-Patell d2002l) studied 
Neon abund a nces o f H n regions in M33, including BCLMP 302. 
iRubin etal] d2008l) used Spitzer-IRS to map the emission 
lines of [Siv] 10.51, H(7-6) 12.37, [Nen] 12.8, [Nem] 15.56 
and [Sm] 18.71 //m in 25 Hn regions in M33, including 
BCLM P 302. ISO/LWS [C n] unresolved spectra at ~ 70" res- 
olution dGrv et al. 1120031) are available for this region from the 
archive. Here, we present PACS maps of [Cn] and [Oi] 63 /mi 
at a resolution of 12", together with a HIFI [Cn] spectrum at 
2 km s _1 velocity resolution. We compare these data with (i) the 
Ha emission tracing the ionized gas, (ii) dust continuum images 
at mid- and far-infrared wavelengths observed with Spitzer and 
Herschel, tracing the dust heated by newly formed stars and the 
diffuse interstellar radiation field, and (iii) CO and H i emission 
tracing the neutral molecular and atomic gas. 

The rest of the paper is organized as follows: Sec. 2 presents 
details of our observations and ancillary data, Sec. 3 states the 
basic results of the spectroscopic observations and a qualitative 
and a quantitative comparison of the [Cn] and [Oi] emission 
with all other available tracers and their correlation. Sec. 4 stud- 
ies the role of [C n] as an indicator of the star formation rate 
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Fig. 3. PACS spectra of [Cn] (158 fim), [Oi] (63 fim), [Nn] 
(122 fim) and [Nm] 57 fim at the Ha peak position of 
BCLMP 302. The LSR velocity is given. All lines are unre- 
solved, i.e. line profiles only reflect the instrumental profiles. 
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Fig. 4. All four spectra of [Cn], Hi and the (2-1) transition of 
CO and I3 CO at the Ha peak position of the H n region BCLMP 
302. All four spectra are at ~ 12" resolution allowing for a 
detailed comparison. The vertical lines denote the velocities - 
249.7, -252.8, -254. lkms' 1 corresponding to H i, CO, and [C n] 
emission respectively. 



2. Observations 

2.1. Herschel: PACS Mapping 

A region extending over 2'x2' around the Hn region BCLMP 
302 in the northern arm of M 33 was observed with the 5x5 pixel 
integral field unit (IFU) of the PACS Spectrometer using the 
wavelength switching (WS) mode in combination with observa- 
tions of an emission free off source position outside of the galaxy 
at PvA/Dec (J2000) = 22.587173 0.6404°. The field of v iew of the 
IFU is 47"x47" with 9'.'4 pixels dPoglitsch et alj2010l) . We used 
the 1st and 3rd order gratings to observe the [Cn] 157.7 /urn, [Oi] 
63. 18 /an, [Oi] 145.52 /an, [Nn] 121.9 /an, [Nm] 57.3 /an, and 
[N n] 205. 18 fim lines, with the shortest possible observing time 
(1 line repetition, 1 cycle), and with a reference position at RA = 
01 h 30 m 20 s .9, Dec = 30°38'25'.'4 (J2000). The reference position 
was selected based on Hi and 100//m ISSAQmaps. 

For all the lines a 3x3 raster was observed on a 40" grid 
with the IFU centered at R.A. = 01 A 34 m 09.9 Dec=+30°47'18'/6 
(J2000) with position angle, PA. = 22.5°. The resulting footprint 
is shown in Fig. [2] The FWHM beam size of the PACS spec- 
trometer is 9.2" near 63 /an and 11.2" near 158/an (E.Sturm, 
priv. comm.). The lines are unresolved, as the spectral resolu- 
tion of PACS is larger than 90km s~' for all lines. The obser- 
vations were performed on January 7, 2010 and the total ob- 
serving time was 1.1 hours for all the 6 lines. T he PACS spectr a 
were reduced using HIPE version 3.0 CIB 1452 (lOtt et al.l2010h . 
The WS data reduction pipeline was custom-made by the NASA 
Herschel Science Center (NHSC) helpdesk. The data were ex- 
ported to FITS cubes which were later analyzed using internally 
developed IDL routines to extract the line intensity maps. 

Using PACS, the [Cn]158/an, [Oi] 63 /an, [Nn] 122 /an, 
and the [N m] 57 fim lines have been detected . The lines of [O i] 
145 fim and the [Nn] 205 fim were not detected. The peak and 
lcr noise limits of the intensities of the [C n], [O i], [N n] 122 fim 
and [N in] 57 fim spectra, at the Ha peak position, are presented 
in TableQ] Fig.|3]show the observed PACS spectra at the position 
of the Ha peak position. The peak integrated intensities were de- 
rived by first fitting and subtracting a polynomial baseline of 2nd 
order and then fitting a Gaussian. Due to unequal coverage of 
different positions, as seen from the grid of observed positions 
shown in Fig. [2] the rms achieved is not uniform. It varies by 
about a factor of 3 over the entire map. 



Table 1. Peak, signal to noise ratio (SNR), and sigma values of 
the integrated line intensities at the Ha peak position for the lines 
detected with PACS all at original resolution. 



Line 




Peak (SNR) 


lcr 






erg s" 1 crrr 2 sr~' 


erg s -1 crrr 2 sr -1 


[Cn] 


158/an 


1.18X10- 4 (67) 


1.72x10-° 


[Or] 


63/jm 


7.20xl(T 5 (6) 


1.20xl(r 5 


[Nn] 


122 fim 


9.95x10-° (7) 


1.42x10"° 


[Nm 


57 fim 


1.30xl(T 5 (7) 


1.80x10-° 



(SFR) and Sec. 5 analyzes the energy balance in the mapped re- 
gion. Sec. 6 presents a detailed analysis of the emission from the 
Ha peak position in BCLMP 302 in terms of models of PDRs. 
In Sec. 7 we summarize and discuss the major findings of the 
paper. 



2.2. Herschel: HIFI spectrum at the Ha peak 

Using HIFI we have observed a single spectrum at the peak po- 
sition (R.A. = 01 h 34 m 06*79 Dec = 30°47'23'.'l (J2000)) of the 
Ha emission from BCLMP 302 (Fig. [4]). The HIFI spectrum was 



1 IRAS Sky Survey Atlas 
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taken on 01 August 2010 during one hour of observing time us- 
ing the load chop mode with the same reference position as was 
used for the PACS observation. The frequency of the [C n] line 
is 1900536.9 MHz, known to wi thin an uncertainty of 1.3 MHz 
(0.2 km s _1 ) dCooksv et alJI 1986b . The blue shifted line required 
to tune the local oscillator to 1899.268 GHz, about the highest 
frequency accessible to HIFI. The [C n] spectra were recorded 
using the wide band acousto optical spectrometer, covering a 
bandwidth of 2.4 GHz for each polarization with a spectral res- 
olution of 1 MHz. We calculated the noise-weighted averaged 
spectrum, combining both polarizations. A fringe fitting tool 
available within HIPE was used to subtract standing waves, sub- 
sequently the data were exported to CLASS for further analysis. 
Next, a linear baseline was subtracted and the spectrum was re- 
binned to a velocity resolution of 0.63 km s _1 (Fig. [4]). We scaled 
the resulting data to the main beam scale using a beam efficiency 
of 69%, using the Ruze formula with the beam efficiency for 
a perfect prima ry mirror ?7 m b .o = 0.76 and a surface accuracy 
of cr = 3.8 /mi (Olberg 20101) . The measured peak temperature 
from a Gaussian fit is 1.14K and the rms (l<x limit) is HOmK 
at 0.63 kms -1 resolution, which is consistent within 10% of the 
rms predicted by HSPOT. The half power beam width (HPBW) 
is 12.2". 



2.3. I RAM 30m CO observations 

For comparison with the [C n] spectrum we have observed spec- 
tra of the (2-1) and (1-0) transitions of 12 CO and I3 CO, at the 
position of the Ha- peak, using the IRAM 30 m telescope on 21 
August 2010. These observations used the backend VESPA. The 
spectra were smoothed to a velocity resolution of 1 km s _1 for 
all the CO lines. The forward and beam efficiencies are 95% and 
80% respectively for the (1-0) transition. The same quantities 
are 90% and 58% for the (2-1) transition. The half -power beam 
widths for the (1-0) and (2-1) transitions are 22" and 12", re- 
spectively. 

2.4. Complementary data for comparison 

We compare the Herschel data of the BCLMP 302 region, with 
the Hi VLA and CO(2-l) HERA /30m map at 12" resolution 
presented by iGratier et alJ (l2010h . We refer to the latter paper 
for a presentation of the noise properties. 

We als o use the Sp it zer M IPS 24 /mi map pre- 
sented by iTabatabaei et aj.l (120071) . the KPNO Ha map 
dHoopes & Walterbosl |2000T) . and the HerM33e s PACS and 
SPIRE maps at 100, 160 , 250, 350 and 500 //m dKramer et alJ 
120 lOt IVerlev et all 1201 Ot iBoquien et all l2010ah . The angular 
resolution of the 100 and 160 //m PACS maps are ~ 6" and 12". 
The rms noise levels of the PACS maps are 2.6 mJy pix 2 at 
100 /mi and 6.9 mJy pix~' at 160 /mi where the pixel sizes are 
3"2 and 674 respectively. 

[C n] observations at two positions within our mapped re- 
gion were extracted from the ISO/LWS archival data. The two 
positions are at RA= 01 h 34 m 07 s Dec=30°46'55" (J2000) and 
RA=01 h 34 m 09 s DEC=30°47'41" (J2000). 



3. Results 

3. 1 . PACS Far-Infrared Spectroscopy 

Within the 2'x2' region mapped with PACS (Fig. 13, we have de- 
tected extended [C n] emission from (a) the northern spiral arm 
traced e.g. by the 100 pm emission, with the strongest emission 



arising towards the H n region BCLMP 302 and (b) from the dif- 
fuse regions to the south-east and north-west. Comparison of the 
[C n] PACS intensities with the ISO/LWS [C n] data at the two 
positions shows an agreement of better than 11% at both posi- 
tions. For the comparison of the intensities at the LWS positions, 
we first convolved the PACS [C n] map to the angular resolution 
of the LWS data. The full ISO/LWS [C n] data set along the ma- 
jor axis of M33 will be published in a separate paper by Abreu 
et al. (in prep.). 

Some of the PACS [O i] 63 /mi spectra displayed baseline 
problems, attributed to the now decommissioned wavelength 
switching mode. Spectra taken along the [C n] emitting part of 
the spiral arm extending north-east to south-west showed prob- 
lems and have been blanked. Both, the [Nn] 122/im and the 
[N m] 57 fim lines were detected at only a few positions within 
the mapped region. The maps of [O i] 63 /mi and [C n] 158 /mi 
emission towards the H n region look very similar and both peak 
at RA = 01 h 34 m 06 s .3, Dec = 30°47'25730. In addition, the [Oi] 
63 /mi map shows a secondary peak towards the south-west, 
to the south of the [Cn] ridge, at RA=01 h 34 m 04 s .364 Dec = 
30°46'26'/55 (J2000), which is not found in the [Cn] map. The 
second [O i] peak lies between the two ridges detected in CO(2- 
1) and coincides with an Hi peak. This suggests that the [Oi] 
emission at this peak position arises in very dense atomic gas. 

Overlays of the [Cn] map with maps of Ha, Hi, CO(2-l) 
and dust continuum in the MIR and FIR (MIPS 24 /mi, PACS 
100 & 160 /mi) are shown in Fig. [5] The dust continuum maps 
correlate well with the [Cn] map. They peak towards the Hn 
region and show the spiral arm extending from the H n region 
in south-western direction. In contrast, CO emission shows a 
clumpy structure wrapping around the H n region towards the 
east. CO emission shows the spiral arm seen in the continuum, 
but its peaks are shifted towards the south. The Hi emission 
shows a completely different morphology, peaking towards the 
north and south of the H n region and showing a clumpy fila- 
ment running towards the west. Further below, we will discuss 
the correlations in more detail. 

Fig. |6]shows overlays of the [O i] 63 /mi map at a resolution 
of 12" with Hi and CO(2-l). Towards the south and south-west 
of the H ii region, the [O i] 63 /mi emission matches the H i emis- 
sion well. This is surprising given the high excitation require- 
ments for the [O i] line. The secondary [O i] 63 fim peak towards 
the south-west, is also traced by H i. It lies in between two ridges 
of CO emission, the arm running north-east to south-west, and a 
second ridge of emission running in north-south direction. The 
northern part of this second ridge shows an interesting layering 
of emission: both [Oi] and Hi are slightly shifted towards the 
east relative to this CO ridge. However, there is much less corre- 
spondence between [O i] and H i emission towards the east and 
north of the H n region. 

3.2. HIFI spectroscopy of the Ha region 

Fig. [4] shows the velocity resolved 158 /mi [Cn] spectrum ob- 
served with HIFI at the Ha peak position. In addition, we show 
the spectra of Hi and the J=2-l transitions of CO and 13 CO. 
HIFI and PACS integrated intensities agree very well. The [C n] 
integrated intensity of 15.6 K km s , observed with HIFI corre- 
sponds to an intensity of l.lOxlO -4 erg s cirT 2 sr~' and this 
matches extremely well with the [C n] intensity observed with 
PACS, at the nearest PACS position, which is offset by only 3". 

All spectra are at a common resolution of ~ 12", allow- 
ing for a detailed comparison. All spectra show a Gaussian line 
shape. However, the line widths are strikingly different between 
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the atomic material traced by H i, the [C n] line, and the molec- 
ular gas. The Hi spectrum shows a FWHM of 16.5 km s , that 
of [Cn] is 13.3 km s _1 , while CO 2-1 shows a width of only 
8 km s _1 (Table |2]i. We take this as an indication that the H i disk 
along the line of sight is much thicker than the molecular disk 
while the material traced by [C n] appears to lie in between. In 
addition, we find that the lines are not centered at the same ve- 
locity. The HI line is shifted by +4.4 km s _1 relative to [Cn], 
while the CO lines are shifted by +1.6 km s relative to [Cn], 
confirming that all three tracers trace different components of the 
ISM. The shifts are significant, as the error of the Doppler cor- 
rections are much smaller, for HIFI (D.Teyssier, priv. comm.) as 
for the other data. 

Based on lower spectral resolution (6 km s -1 ) optical spec- 
troscopy of HB and [Om] emission lines, the velocity of the 
ionized gas is deduced to be shifted by -16.9 km s with re- 
spect to the [C n] line, with a velocity dispersion of 1 1 km s _1 
dWillner & Nelson-Pate1l2002t IZaritskv et al.lll989l) . The optical 
spectroscopic data is however at much higher angular resolution 
(2-4"). A map of the emission lines originating solely from the 
ionized medium, which would allow to smooth these data to 12" 
resolution for direct comparison with the other data, is not yet 
available. In Section 6.2, we will use PDR models to show that 
about 20-30% of the observed [C n] stems from the ionized gas 
of the BCLMP 302 H n region. 

In summary, we find a "layering" of the line-of-sight veloci- 
ties tracing the different ISM components. The velocity increases 
from Ha to [C n] to CO to Hi. 

Table 2. Parameters derived from Gaussian fits to spectra ob- 
served with HIFI, VLA, and IRAM 30m, at the Ha peak position 
(Figs.|2]&|4]i. #b indicates the half power beamwidth. 



Line 


8 b 

tf 


J rdv 
Kkms-' 


^ 7 cen 

km s _I 


Av 
km s 


[Cn] 


11 


15.6+0.5 


-254.1+0.1 


12.9+0.9 


Hi 


11 


1446.7±44.1 


-249.7±0.2 


16.5+0.6 


CO (1-0) 


22 


3.2±0.2 


-252.3±0.2 


7.5+0.5 


CO (2-1) 


11 


5.3+0.1 


-252.8+0.1 


7.9+0.3 


13 CO(1-0) 


22 


0.28±0.02 


-252.9+0.2 


5.8+0.5 


13 CO(2-l) 


11 


0.69+0.06 


-252.8+0.3 


6.2+0.5 



3.3. Detailed comparison of [Cu] emission with other tracers 
Table 3. Correlation coefficients for the scatter plots in Fig. [7] 



Tracers 




Correlation Coefficient 






Entire 


Region 


Region 


Region 




Map 


A 


B 


C 


[Cn]-Ha 


0.60 


0.93 


0.66 


0.30 


[Cn]-CO(2-l) 


0.41 


0.40 


0.47 




[Cn]-Hi 


0.22 


< 0.1 


<0.10 


0.14 


[CiiHOi] 




0.77 






[Cn]-F 24 


0.61 


0.86 


0.66 


0.41 


[Cn]-F 100 


0.59 


0.85 


0.56 


0.43 



For a more quantitative estimate of the correspondence be- 
tween the different tracers in which the spiral arm has been 



mapped we have created scatterplots of intensities of tracers like 
Ha, 12 CO, Hi, [Oi] 63 jum, MIPS 24 /an and PACS 100 /on as 
a function of the [Cn] intensities (Fig.|7]i. We have used inten- 
sities from all the maps which are smoothed to a resolution of 
12" and gridded on a 12" grid. In order to identify any apparent 
trends in the emission we have defined three sub-regions within 
the mapped region: Region A corresponds to the Hn region 
BCLMP 302 itself, Region B corresponds to the south-western 
more quiescent part of the spiral arm, traced by e.g. the 100 /an 
continuum emission, and centered on a peak of CO emission, 
and Region C lies outside of the prominent CO arm. These three 
sub-regions are marked by black rectangles in Figure[5] A sec- 
ond H ii region along the arm, lies just outside and to the north of 
the box defining the quiescent arm region. For the remainder of 
the paper we always analyze and compared the results for these 
three sub-regions separately. 

In the log-log scatter plots of Figure|7] we find that the Ha, 
[Oi], and continuum emission at 24pm and 100 /on show pro- 
nounced linear correlations with the [C n] emission within the 
H n region. In the region C, the intensities of the Ha emission 
and the continuum emission at 24 and 100 /mi remain almost 
constant. The CO(2-l) intensity in the Hn region (A) is only 
poorly correlated with the [Cn] emission. In regions B and C 
the CO(2-l) intensity shows no correlation with the [C n] inten- 
sity and has a large scatter. Hi does not show any correlation 
with [Cn]. 

A more quantitative analysis of the correlation between the 
different tracers and [C n], is obtained by calculating the correla- 
tion coefficient (r) (Tabled. For the entire region Ha and 24 fim 
and lOOjt/m intensities show around 60% correlation with the 
[Cn] intensities. For the region A Ha, [Oi], 24 //m, and lOOyt/m 
are well correlated (r > 0.75) with [Cn]. The Ha emission is 
strongly correlated with the [Cn] intensity also on the south- 
western arm position. The [Oi]/[Cn] intensity ratio measured 
primarily in the H n region varies between 0. 1-0.4 and this vari- 
ation is significantly larger than the estimated uncertainties. 



4. [C n] as a tracer of star formation 

In Fig. [8j we plot the star formation rate (SFR), estimated 
from the 24 fim MIPS data and the KPNO Ha data, as a func- 
tion of the [Cn] intensity. Positions within the selected re- 
gions A, B, and C are marked using different symbols and 
colours similar to Fig. [7] The SFR has been calculated from 
SFR = [L(Ho) + 0.031L(24)] x 5.35 x 10~ 35 in M yr" 1 , 
where L(Ha) is the Ha luminosity in Watt a nd L(24) is defined 
a s vL v at 24 /an in Watt (ICalzetti e t al. 2007). We have assumed 
a lKroupal ( 120011) initial mass function with a constant SFR over 
100 Myr. Here, we study the correlation on scales of 12" corre- 
sponding to 50 pc. On these small scales, we may start to see a 
break-up of any tight correlations between the various tracers of 
the SFR. 

Viewing at all pixels we find a steepening of the slope in 
the log-log plots from regions where both the SFR and [C n] 
are weak, to regions where both are strong. Towards the Hn re- 
gion (A), we find an almost linear relation between log(SFR) 
and log([Cn]), with a good correlation, r = 0.90 and the fitted 
slope is 1.48+0.26, i.e. the SFR goes as [Cn] to the power 1.48. 
Region B shows a correlation coefficient of 0.64 and region C 
shows no correlation. 
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Fig. 8. Correlation of [C n] intensities with SFR (derived from 
24/im + Ha). Red triangles correspond to positions in region 
A region, the green squares represent positions in region B and 
the cyan filled circles correspond to positions in region C. Each 
marker corresponds to one position on a 12" grid. Details about 
the regions and points are identical to Fig. [7] The black straight 
line marks the fits to the region A. The errorbar corresponds to 
an uncertainty of 20%. 
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Fig. 9. Correlation plot of [C n] and Total Infrared (TIR) inten- 
sities at a 12" resolution and on a 12" grid. Markers and the 
straight line are the same as in Fig. [7] 



5. Energy Balance in the spiral arm 

iBoquien et al.l d2010bl) have obtained a linear fit to the total 
infrared (TIR) luminosity as a function of the luminosity in 
the PACS 160 pm band for the entire M33 galaxy. The TIR 
is the total infrared intensi ty, integrated betw een 1 fim and 
1 mm. It is about a factor 2 dRubin et al.ll2009l) larger than the 
FIR continuum which is integrated between 42.7 and 122 urn. 
IBoquien et all (1201 Obi) derived the TIR from fits of iDraine et"ai1 
(120071) models to the MIPS, PACS, and SPIRE FIR data. And 
they find a tight linear relation between the two quantities: 
logL TIR -ax logLieo + b with a = 1.013+0.008 and b= 
0.429+0.097. We have used this relation to derive the TIR in- 



* 



TIR (erg/s/cm /sr) 



Fig. 10. Plot of [C n]/TIR and ([C n]+ [O i])/TIR (blue asterisks) 
in the H n region as a function of the total infrared (TIR) inten- 
sities at a resolution of 12" on a 6" grid. The error bars denote 
30% errors. Markers same as in Fig. [7] 



tensity at each position of the mapped region, at a resolution of 
12". To check that this is a valid approach, we independently 
derived the TIR at individual positions by fitting a greybody to 
the MIPS, PACS, and SPIRE data, smoothed to a common reso- 
lution of 40". The resulting TIR agrees within 10% with the TIR 
derived from only the 160 fim band. 

Figure [9] shows a scatterplot between TIR and [Cn]. 
Positions corresponding to the selected sub-regions are indicated 
using different markers. We find that [C n] and TIR are tightly 
correlated with a correlation coefficient of 0.87 only in the Hn 
region (A) with the fitted slope being 0.64+0.14. The correlation 
coefficient in region B is 0.52 and and C shows no correlation. 

Incident FUV photons with energies high enough to eject 
electrons from dust grains (hv > 6 eV) heat the gas via 
these photoelectrons, with a typical efficiency of 0.1-1% 
dHollenbach & Tielensl 19971) . Efficiency is defined as the energy 
input to the gas divided by the total energy of the FUV pho- 
tons absorbed by dust gr ains. Based on ISO/LW S observations 
of a sample of galaxies, iMalhotra et al.l (1200 ll) found that (a) 
more than 60% of the galaxies show L[ch]/Lfir > 0.2% and (b) 
L[Cii] /Lhr decreases with warmer FIR colors and increasing star 
formation activity, indicated by higher Lfir/^b ratios, where Lb 
is the luminosity in the B band. We have calculated the ratio of 
the intensities of [C n]/TIR, as a proxy for the heating efficiency, 
and plotted it against with [C n] intensities (Fig. [Tot. The heating 
efficiency vary by more than one order of magnitude within the 
2' x 2' region, between 0.07 and 1.5%. Considering an uncer- 
tainty of 20% in both the measured 160/im intensities and the 
[C n] intensities, for the relation between the 160 fim luminosity 
and the TIR luminosity mentioned earlier, we estimate the un- 
certainty in the [C n]/TIR ratio to be ~ 30%. The observed vari- 
ation in the [C n]/TIR intensity ratios is significantly larger than 
the uncertainty we estimate from the [Cn] and 160/um intensi- 
ties. For the H n region (A), we find heating efficiencies between 
0.2 and 1.0%. Regions B and C show efficiencies between 0.15- 
0.4% and 0.4-1.2%. Within the Hn region (A), the total heat- 
ing efficiency, including the [Oi] 63 pm line, ([Cn]+[Oi])/TIR, 
lies between 0.3-1.2% (Fig.lTOll. Outside the Hn region, reliable 
[O i] data is largely missing. 
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Table 5. Properties of the H n region BCLMP 302 



Excitation parameter u 180 pc cm 1 IK74 

Radius r 39 pc IK74 

Mass 10 5 M IK74 

rms electron density 6.2 cirT 3 IK74 

Electron density n e 100 cm -3 R08 

Ionization parameter U -3.3 this paper 

Effective temperature r c) f 38,000 K this paper 



Ho- luminosity 2.2 10 38 erg s -1 KPNO map 

Notes. As explained in the text, the value given for the electron density 
(n e )is assumed to hold for BCLMP 302. 

Refer ences. IK74 llsrael & van der Kruitl dl974l) . R08 iRubin et alj 
d2008h 



6. Modeling the PDR emission at the Ha peak of 
BCLMP 302 

6.1. Estimate of FUV from TIR 

One of the key parameters of PDRs is the FUV radiation 
field which heats the PDR. From energy considerations, the 
total infrared cooling emission is a measure of the irradiat- 
ing FUV photons of the embedded OB stars. We estimate 
the FUV flux Go (6eV< hv <13.6eV) impinging onto the 
cloud surfaces from the emergent total infrared intensities: Go 
— An /fir - An 0.5 / tir with Go in units of the Habing field 
1.6 10~ 3 ergs -1 cm -2 dHabind fl968h and the intensity of the 
far-infrared continuum between 42.5 //m and 122.5 /mi, /fir, 
in units of erg s -1 cirT 2 sr~'. Here, we assume that heating 

by photons with hv < 6eV contributes a factor of ~ 2 

dTielens & Hollenbach|[l985l) and that the bolometric dust con- 
tinuum intensity /tir is a factor of ~ 2 larger than /fir dDale et al.l 
l200l . 

At the Ha peak position, the TIR intensity of 
1.18x10 2 erg s -1 crrT 2 sr -1 translates into a FUV field of 
Go = 46 in Habing units. Outside of the Ha peak, the FUV 
field, estimated from the TIR, drops by more than one order of 
magnitude (cf. Fig. [9). 

The FUV radiation lea king out of the clouds is measured 
using the GALEX UV data (iMartin et alJl2005t iGif de Paz et al.l 
120071) to be G = 24 for a 12" aperture. 

We thus find, that 66% of the FUV photons emitted by the 
OB stars of the H 11 region are absorbed and re-radiated by the 
dust, and 34% are leaking out of the cloud, at the Ha position. 
This is consistent with the FUV extinction de rived from the Ha 
and 24 //m fluxe s dRelafio & Kennicuttll2009l) and the FUV/Ha 
reddening curve dCalzettill200iF ~ 

It is interesting to note that Hn regions observed in the 
LMC and other spi ral galaxies by lOev & Kennicutj d 19971) and 
iRelano et alj (Hp02) show that typically around 50% of the ion- 
izing stellar radiation escape the H 11 regions, roughly similar to 
the fraction of 36% we find in BCLMP 302/M33. 

6.2. [Cn] emission from the ionized gas 

In Table we present properties of the H n region BCLMP 302 
compiled from the literature and also calculated thereof. These 
properties were used to estimate the fraction of [C 11] emission 
contributed by th is H n region using the model calculations of 
lAbel et al.l d2005l) and lAbell d2006l) . These authors have esti- 
mated the [C n] emission for a wide range of physical conditions 



in H 11 regions, varying their electron density n e , ionization pa- 
rameter U, and effective temperature r e jf . 

The [C ill PACS map is centere d on the bright Hn region 
BCLMP302 dBoulesteix et alj|1974h . wh ich corresponds to the 
H n re gion no. 53 in the M33 catalog of llsrael & van der Kruitl 
(1 1974 IK74). Using the measured radio fl ux and the formul a 
M (pccirr 2 ) = 13.5(5 /f.u.) 1/3 (D/kpc) 2/3 dlsrael et al.l [19731 
IK74 estimate the excitation parameter u = 180 pc cirT 2 for a 
distance of 720 kpc. IK74 also derive an rms electron density 
of 6.2cm 3 , and estimate the radius r of the Hn region to be 
38.5 pc. We did not correct these results for the now better known 
distance, as it does not affe ct our conclusio ns. For an electron 
temperature T e of 10,000 K. iPanagiald 19731) expresses the exci- 
tation parameter as «(pc cirT 2 ) = 2.2xl0~ 19 In(L)(J3 - /JiT 1 ]'^, 
where ifi-fii) is the recombination rate to the excited levels of 
hydrogen in units of cm 3 s -1 . Thus we get for the total flux of 
ionizing photons N(L) = 4.6 10 43 m 3 s~'. Using the above val- 
ues of the parameters, the dimensionless ionization parameter U 
was derived via t/ = MX)/(c47rr 2 n H+ )dl vans & Dopita|[l98l 
lMorissetll2004l) . 

The rms electron density derived by IK74, provides only a 
lower limit to the true electron density and observations of e.g. 
the [S n] doubl et at 6754 A are m issing which would allow a 
direct estimate. IRubin etaO d2008l) observed 25 H n regions in 
M33 using Spitzer-IRS, and concluded that their electron den- 
sities are ~ 100cm 3 . Since there is no apparent reason for 
the n e in BCLMP 302 to be significantly different, we assume 
n e = 100 cirT 3 . Hence, with «h+ = n e , we estimate log U = -3.3 
(Table©. 

The observed ratio of double and single ionized Nitrogen, 
[Nm] 57/mi/[Nii] 122 /mi, at the Ha position (Table[TJ, indi- 
cates an effective temperat ure of the ionizing stars of about 
38.000K dRubin et al.lll994l) . 

From the model calculations of lAbell (120061 Fig. 5), we es- 
timate the fraction of [C 11] emission from the BCLMP 302 H 11 
region. The [C 11] fraction increases with dropping electron den- 
sity, but is only weakly dependent on the ionization parameter 
and the stellar continuum model. Depending on the model, the 
resulting fraction varies slightly, ~ 20-30%. The fraction would 
drop to about 10%, for an electron density of 10 3 cirT 3 . 

Observations of the [N n] 205 fim line, in addition to the 
122 /mi line would allow to estimate more accurately the elec- 
tron density of the H 11 region. The model calculations and ob- 
servations compiled by lAbell d2006l) show that the [C 11] intensi- 
ties stemming from H 11 regions and the intensities of the [N 11] 
205 fim line are tightly correlated. lAbell d2006l) find Log [/£!] = 
0.937Log[/™ 1 ] + 0.689 (erg cirT 2 s _1 ). Assuming that 30% 
of [Cn] emission stems from the BCLMP 302 Hn region at 
the Ha peak position, we estimate a [N 11] 205 fim intensity of 
~ 3.6 x 10~ 6 erg s~' cirT 2 sr~' (cf. TableQ}, corresponding to 
3 x 10~ 18 WirT 2 . The expected [Nn] 205 /mi intensity is a fac- 
tor of 3 below the estimated rms of our PACS observations and 
hence is consistent with the non-detection. Using the detected 
[N 11] 122 fim line, and the predicted [N 11] 205 fim intensity from 
above, the [N 11] 122 //m/[N 11 ] 205 um ratio is 2.8 and this corre- 
sponds to an n e of 10 2 cirT 3 dAbel et al. 2005, Fig. 22), which is 
consistent with the n e assumed by us for the above calculations. 

6.3. Comparison with PDR models 

Here we investigate whether the FIR and millimeter lines, to- 
gether with the TIR, observed towards the Ha peak, can be con- 
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Table 6. Input to the PDR model for the HIFI position 



n [cm ] 



Fig. 11. Comparison of line intensities and intensity ratios at 
the position o f the Ha peak with pl ane-parallel constant density 
PDR models (iKaufman et al.llT999b . Grey-scales show the esti- 
mated reduced^ 2 . The horizontal line shows the FUV estimated 
from the total FIR intensity. The contours correspond to the in- 
tensities/ratios of different spectral lines (as shown in the labels). 
The [C n] intensity corresponds only to the estimated contribu- 
tion (70% of the total) from the PDRs. 
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CII/CO(1-0) 

Fig. 12. Star formation rate (SFR) versus the [C n]/CO(l-0) in- 
tensity ratio. Markers are the same as in Figure|7] 



sistently explained in terms of emission from the PDRs at the 
surfaces of the molecular clouds. We use the line and total in- 
frared continuum intensities at a common resolution of 12". 

At the position of the Ha peak we observe the follow- 
ing intensity ratios on the erg scale (cf. Tabled [Oi]/[Cn] 
= 0.35, [Cn]/CO(l-0) = l.lxlO 4 , [Cn]/CO(2-l) = 1263, 
([Cii]+[Oi])/TIR = 9.6x10 3 and CO(2-1)/C O(1-0) = 8.8 
Comp aring all ratios with the PDR model of IKaufman et al] 
(fl999h (Fig.HB, we find a best fitting solution near a FUV field 
of Go=32 in units of the Habing field and a density of 320 cm -3 . 
The fitted FUV field agrees rather well with the FUV field of 
Go = 46 (in Habing units) estimated from the TIR continuum. 
The corresponding reduced^ 2 was estimated assuming an error 



Tracer 


Intensity 


[Cn]' 


8.4 lO" 5 


[Oi] 


3.0 10- 5 


co(i-or 


7.5 10-* 


CO (2-1) 


6.7 10~ 8 


TIR 


1.18 10~ 2 



Notes. All intensities are in units of erg sec~'cnT 2 sr~'. * The observed 
[C n] intensities was multiplied by 0.7, assuming that 30% of the emis- 
sion stems from the H n region. ** The CO 1-0 data at 22" resolution 
were multiplied by the beamfilling factor 1.5, estimated by smoothing 
the 2-1 data from 11" to 22" resolution. 



of 30% on the observed intensity ratios. The reduced x is de- 
fined as x 1 = !/(« - 2)(Z (I? s ~ I"' oM ) 2 /(T 2 ), where (« - 2) is 

the number of degrees of freedom with n being the number of 
observed quantities used for the fitting and / denotes either inte- 
grated intensities or ratios of integrated intensities. We find the 
minimum value of x 2 to be 4.4, indicating that the fit is not sat- 
isfactory. Indeed, the fitted density seems rather low, given the 
high critical densities of the [C n] and in particular also of the 
[Oi] line. On the other hand, the observed [Oi]/[Cii] ratio, is 
consistent with the low density solution. The ratios with the CO 
lines, deviate from this solution. For instance, the CO 2-1/1-0 
ratio together with the derived FUV, indicates a higher density 
of about 10 4 cirT 3 . Using line ratios to compare with the PDR 
model, allows to ignore beam filling effects, to first order. Note 
that the absolute [C n] intensity, reduced by the fraction of 30% 
stemming from the H n region, agrees well with the best fitting 
solution, indicating a [C n] beam filling factor of about 1 . 

The rather poor fit of the intensity ratios towards the Ha peak 
shows the short comings of a plane-parallel single density PDR 
model . First tests using KOSMA-r PDR models iRollig et all 
(120061) of spherical clumps with density gradients reconcile bet- 
ter to the observations. This indicates strong density gradients 
along the line of sight. In a second paper, we shall include new 
HIFI [C n] data along two cuts through the BCLMP 302 region 
and explore detailed PDR models, which include the effects of 
geometry and sub-solar metallicity. 



7. Discussion 

Mapping observations of the northern inner arm of M 33 at an 
unprecedented spatial resolution of 50 pc have revealed details 
of the distribution of the various components of the interstel- 
lar medium and their contribution to the [C n] emission. [C n] at 
158 fim is one of the major cooling lines of the interstellar gas. 
Thus irrespective of whether hydrogen is atomic or molecular, 
the [C n] line emission is expected to be strong wherever there 
is warm and photodissociated gas. We have identified emission 
towards the H n region as well as from the spiral arm seen in the 
continuum, and from a region outside of the well-defined spi- 
ral arm. [C n] emission is strongly correlated with the Ha and 
dust continuum emission, while there is little correlation with 
CO, and even le ss with H i. This su ggests that the cold neutral 
medium (CNM; IWolfire et al.ll 19951) does not contribute signifi- 
cantly to the JC n] em ission in the BCLMP 302 region. Recently, 
lLanger et al.l (1201 01) found a similar poor correlation between 
[C n] and H i emission in a sample of 29 diffuse clouds using 
HIFI. The lack of correlation between [C n] and CO found in 
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BCLMP 302, may indicate that significant parts of the molecu- 
lar gas are not traced by CO because it is photo-dissociated in the 
low-metallicity environment of M33. This interpre tation is con- 
sistent with both theoretical models developed by iBolatto et al.l 
(fl999h and recent observational studies of diffuse clouds in the 
Milky Way by Langer et al.l d2010l) . and of dwarf galaxies by 
IMadden et aTT(l201 lb . 



Comparison of the first velocity-resolved [C n] spectrum of 
M33 (at the Ha peak of BCLMP 302) with CO line profiles show 
that the [C n] profile is much broader, by a factor of ~ 1 .6, and 
slightly shifted in velocity, by ~ 1 .6 km s _1 . Compared to the H i 
line, at the same angular resolution, the [C n] line is less broad 
by a factor ~ 1.3 and shifted by ~ 4.4 km s . These findings 
indicate that the [C n] line is not completely mixed with the CO 
emitting gas, but rather traces a different more turbulent outer 
layer of gas with slightly different systemic velocities, which is 
associated with the ionized gas. 

Interestingly, recent [C nl HIFI observations of Galacti c star 
forming regions dOssenkopf et alJl2010t iJoblin et"al1l2010h also 
show broadened and slightly shifted [C n] line profiles relative to 
CO. 

The two major cooling lines of PDRs are the [Oi] line at 
63 //m and the [C n] 158 fim line. The intensity ratios [OI]/[CII] 
and ([OI]+[CII]) vs. the TIR continuum, have been used exten- 
sively to estimate the de nsity and FUV field of the emitting re- 
gions. Using ISO/LWS, iHiedon et all (120031) observed the far- 
infrared spectra of the nucleus and six giant H n regions in M 33, 
not including BCLMP 302, but including NGC 604, IC 142, and 
NGC595 shown in Figured] The 70" ISO/LWS beam corre- 
sponds to 285 pc and therefore samples a mixture of the differ- 
ent ISM phases. They find [Oi]/[Cii] line ratios in the range 
0.7 to 1.3, similar to the range of values found in the cente r 
and spiral arm positions of M83 and M51 (iKramer et al 1 l2005h . 
Towards the H n region BCLMP 302 in the northern inner arm 
of M 33, we measure much lower [O i]/[C n] ratios between 0.1- 
0. 4. These ratios l i e towa rds the lower end of the values found 
bv lMalhotra et al.l d2001l) in their ISO/LWS study of the unre- 
solved emission of 60 galaxies who find values between ~ 0.2 
and 2. Similarly low values of down to 0.16 are foun d e.g. in the 
Galactic star forming reg ions DR21 and W3IRS5 dJakob et alj 
20071: IKramer et al.ll2004t) . Comp arison with the PDR models of 



Kaufman et al 



d 19991 Fig. 4) and lRolhget"aT1 d2006l) show that 
the [O i] 63 fim line becomes stronger than the [C n] emission in 
regions of high densities of more than about 10 4 cm -3 . At the Ho- 
peak position in BCLMP 302, we observed a ratio of 0.4, after 
correcting [Cn] for the contribution from the ionized gas. This 
ratio indicates lower densities and a FUV field of less than about 
100 Go (cf. Fig-fTTTi. Still lower ratios, indicate lower impinging 
FUV fields. 

The ratio of [Oi] +[Cn] emission over the FIR continuum, 
is a good measure of the total cooling of the gas relative to the 
cooling of the dust, reflecting the ratio of FUV energy heating 
the gas to the FUV energy heating the grains, and hence the 
grain heatin g efficiency, i.e. th e efficiency of the photo-electric 
(PE) effect dRubin et al.ll2009h . Efficiencies of up to about 5% 
are stil l consistent with FUV heating, i.e. with emis sion from 
PDRs (lBakes&Tielenslll994t iKaufman et allll999l) . The PE 
heating efficiency is a function of FUV field, electron density, 
and temperature. A high FUV field leads to a large fraction of 
ionized dust particles, lowering the efficiency. On the other hand, 
low metallicities naturally lead to increased efficiencies as e.g. 
the PDRs become larger when the dust attenuation is lowered 
dRubin et al.ll2009h . 



iHigdon et al.l d2003l) compared the [O i] + [C n] emission with 
the FIR(LWS) continuum integrated between 43 and 197/iin, 
and obtained ratios of 0.2% to 0.7%, corresponding to 0.4 to 
1 .4% for a rough estimate of the FIR /TIR conversion factor of 2 
dDale et al.ll200ll:lRubin et alj|2009l) . Towards the 2' x 2' region 
presented here, the [C n]/TIR ratio varies by more than a factor 
10 between 0.07 and 1.5%. Extra galactic observations at reso- 
lutions of 1 kpc or more, find efficiencies of only up to ~ 0.3% 
dKramer et al.ll2005t iMalhotra et alj|2001l). T he observations of 
M 31 bv iRodriguez -Fernandez et all d2006h with ISO/LWS at 
~ 300 pc re s olutio n, do also find high efficiencies of up to 2%. 
iRubin et alj d2009l) analyzed BICE [Cn] maps of the LMC at 
225 pc resolution, and find efficiencies of upto 4% in this low 
metallicity environment. The mean values found in the LMC are 
however much lower: 0.8% in the diffuse regions, dropping to 
~ 0.4% in the SF regions. A s imilar variation was fou nd in the 
Galactic plane observations by Nakaga wa et al.l d!998l) . Heating 
efficiencies o bserved in the Milky Way span about two orders 
of magnitude. lHabart et al.l d200ll) observed efficiencies as high 
as 3% in the low- UV irradiated Galactic PDR LI 721, whereas 
IVastel et alj 120011) also using ISO/LWS fluxes of [C n] and [0 1], 
found a very low heating efficiency of 0.01% in W49N. The 
global val ue for the Milky W ay from COBE observations is 
~ 0.1 5% (Wright et al. 19911) . The galactic averages observed 
by e.g. IMalhotra et al.l ( 20011) . are dominated by bright emission 
from the nuclei where the SFR and the FUV fields are large, 
hence lowering the photo-electric heating efficiencies. 

Unresolved observations of external galaxies find a tight cor- 
rel ation between CO an d [C 11] emission. This was initially seen 
by IStacev et al.l d!991l) with the KAO, and subsequently sup- 
ported by ISO/LWS observations. However at spatial scales of 
50 pc we do not detect such a tight correlation in BCLMP 302. 
This lack of correlation between CO and [C 11] emission is al- 
ready se en in the maps of spiral arms in M 3 1 at ~ 300 pc res- 
olution dRodriguez-Fernandez et al]|2006l) . As summarized by 
them, the [Cn]/CO(l-0) ratio varies from ~ 1300 in galactic 
disks to about 6000 in starbursts and to about 23,000 in the 
LMC. In BCLMP 302/M 33, we find that the [C n]/CO(l-0) ratio 
varies between 200 to 6000, with the H 11 region having values 
between 800-5000 (Fig. 1121. The region C shows large values 
of the [Cn]/CO(l-0) intensity ratios as CO is hardly detected. 
Thus, while the [C n]/CO(l-0) intensity ratios are higher for re- 
gions with high SFRs, there is no marked correlation between 
the two quantities at scales below about 300 pc. The tight cor- 
relation between [Cn] and CO emission breaks down at scales 
resolving the spiral arms of galaxies. 
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Fig. 5. Overlay of the Ha, Hi and 12 CO(2-l) emission and the dust continuum at 24 (MIPS), 100 & 160/mi (PACS) emission with 
contours of [Cn] at 158/im. The white and black contours are for intensities between 10-20% (in steps of 10%) and 30-90% (in 
steps of 15%) of the peak [Cn] intensity of 1.18xl0~ 4 erg s" 1 cirT 2 sr _1 . The beam sizes for Ha, Hi, 12 CO(2-l), 24 fjm, PACS 
100 fim, PACS 160 fim and 12 CO(2-l) observations are 1", 12", 12", 6", 6'.'7 and 11"4 respectively. The Hi map is integrated 
between -280 km/s to -130 km s _1 . and the CO(2-l) map is integrated between -270 to -220 km s _1 . Marked in the 24jum image 
(with black rectangles) are also the regions selected for further analysis. 
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Fig. 7. Correlation of [Cn] with Ha, 12 CO(2-l), Hi, [Oi] 63 fim, MIPS 24 /mi and PACS 100 /mi. The black crosses mark all 
positions on a 6" grid. The red triangles correspond to positions in region A region, the green squares represent positions in region 
B and the cyan filled circles correspond to positions in region C. All these regions are marked in the 24 fim map in Fig. [5] Errorbars 
corresponding to a 20% error on the plotted quantities along both axes are shown at one representative point in each panel. For 
[0 1] 63 fim, only positions close to the H 11 region are used. The blue dashed line in the [C n]-[0 1] scatter plot corresponds to equal 
intensities of the two tracers. 

Table 4. Average values and their variation of the different tracers arising from the H 11 region, the spiral arm and the region C. 



Tracer 


Region A 


Region B 


Region C 


[Cn] (erg s _1 crrT 2 sr ) 


(4.51 ±2.30)xl0" 5 


(1.28±0.48)xl0" 5 


(1.19 ± 0.42)xl0~ 5 


Ho- (ergs -1 ) 


(7.1±6.9)xl0 37 


(5.3±1.8)xl0 36 


(7.95+1. 93)xl0 36 


CO (2-1) (Kkms- 1 ) 


1.88+0.98 


2.95+1.40 


0.29+0.14 


Hi((Kkms-') 


1597±358 


1267±325 


620+309 


[Oi] (erg s _I cnT 2 sr -1 ) 


(1.28±0.09xl0- 5 






F24 (JybearrT 1 ) 


101.9+82.8 


32.2±8.6 


11.9+1.50 


F100 (JybearrT 1 ) 


0.66±0.37 


0.26±0.05 


0.12+0.03 
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